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Abstract: We study the prospects for detecting a light boson X with mass mx % 
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100 MeV at a low energy electron-proton collider. We focus on the case where X domi- 
1 nantly decays to e + e~ as motivated by recent "dark force" models. In order to evade direct 

^ and indirect constraints, X must have small couplings to the standard model (ax ^ 10 -8 ) 

and a sufficiently large mass (mx > 10 MeV). By comparing the signal and background 
cross sections for the e~pe + e~ final state, we conclude that dark force detection requires 
an integrated luminosity of around 1 ab" 1 , achievable with a forthcoming JLab proposal. 
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1. Introduction 



While the gravitational evidence for dark matter is overwhelming [1, 2, 3, 4], direct mea- 
surements of the spectrum and properties of dark matter have so far been elusive. However, 
recent astrophysical anomalies — including the WMAP Haze [5, 6], the PAMELA, FERMI, 
and H.E.S.S. e+/e~ excesses [7, 8, 9, 10], and the INTEGRAL 511 keV excess [11, 12]- 
could be evidence for dark matter annihilation, decay, or up-scattering in our galactic halo. 
With these observations, an intriguing paradigm for dark matter has emerged, where TeV- 
scale dark matter interacts with a GeV-scale boson [13, 14, 15, 16]. This new light boson 
X typically has a mass in the range 

2m e < m x < few GeV, (1.1) 

with an 0(1) branching fraction X — > e + e - . 1 

What is the best way to look for light bosons with small couplings to the standard 
model? Indirect constraints from lepton anomalous magnetic moments require the coupling 
of X to leptons to be ax < 10~ 8 [18, 19], much smaller than the electromagnetic coupling 
ciem — 1/137. 2 Therefore, any direct production mode for X faces a large irreducible 
background from an equivalent process where X is replaced by an off-shell photon 7*. A 
number of studies at lepton colliders have concluded that around 1 ab _1 of data is needed 
to see the process e + e~ — ► 7 + X [22, 18, 23, 24, 21, 25, 26]. While such large integrated 
luminosities have been achieved at the i?-factories, it is worthwhile to consider alternative 
experimental setups that might be more easily scaled to multi-ab data sets. 

One standard method to find new particles with small couplings is fixed-target exper- 
iments, either with a high intensity beam on a thin target or a "beam dump" experiment 
on a thick target. As we will see, such experiments already constrain the X parameter 
space [27, 28]. Recent studies in Refs. [21, 29] have concluded that improved fixed-target 
experiments can cover a wide range of masses and coupling for X, especially if X has a 
sufficiently long lifetime to yield a displaced vertex or if X has a decay mode to penetrating 
muons. Even in the case of prompt X boson decay, the luminosity achievable in traditional 
fixed-target experiments approaches 1 ab~ 1 /day, so with good energy resolution and con- 
trol over systematics, the irreducible 7* background could be beaten by statistics, and one 
can simply search for electron pairs that reconstruct a narrow X resonance. However, full 
event reconstruction is impossible in this context, since one cannot measure the spectrum 
of the recoiling nucleus, so traditional fixed-target experiments lack a crucial kinematic 
cross-check that is available in lepton colliders. 

In this paper, we propose searching for an X boson in low-energy electron-proton 
collisions through the process 

e'p^ e~p + X, X^e + e~. (1.2) 

With a high intensity electron beam on a diffuse hydrogen gas target, one combines the 
high statistics of a traditional fixed-target experiment with the full event reconstruction 

x For a recent study of models with even lighter bosons, see Ref. [17]. 

2 There are additional direct constraints on X from rare meson decays [18, 20, 19, 21]. 
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potential of a lepton collider. To our knowledge, this experimental setup was first suggested 
in Ref. [30], motivated by a different dark matter scenario with an invisibly decaying X 
boson [31]. In that context, the recoiling proton spectrum was crucial for discovery. Here, 
we focus on X bosons that decay visibly to e + e~. Like the fixed-target proposals in 
Refs. [21, 29], one is still looking for a narrow X resonance on top of a huge radiative 
QED background, but here the recoiling proton and electron spectrum can be used to 
over-constrain the kinematics. 

As in Ref. [30], we consider an electron beam with energy E e ~ 100 MeV, where 
the scattering is dominated by elastic scattering and associated radiative processes. In 
particular, pion production is kinematically forbidden as are nuclear excitations. Such a 
setup is being actively considered for installation at the Free Electron Laser (FEL) at the 
Thomas Jefferson Lab National Accelerator Facility (JLab), replacing the laser cavity with 
a hydrogen gas target. 3 Where the reach in X parameter space overlaps, this proposal of 
a high intensity beam on a diffuse target is complementary to the proposal in Ref. [29] 
of a diffuse beam on a high density target. We will argue that for the same integrated 
luminosity, the X reach in ep collisions is comparable to e + e~ and e~e~ collisions. Since 
1 ab _1 /month is achievable with the FEL beam on a hydrogen gas target, low-energy ep 
collisions are in principle a cost-effective way to search for the X boson. 

In our study, we will focus on irreducible physics backgrounds and assume idealized 
detectors. While there are important experimental backgrounds, we will assume that these 
can be controlled using, for example, information about recoiling proton and electron. The 
JLab FEL setup is in principle sensitive to: 

mx < 2m e : X — > 77, invisible 
2m e < m x < 100 MeV : X -> e+e", 77, invisible (1.3) 

In models with rich dark sectors, one can even imagine multi-body X decay modes or 
more than one X field [33, 24]. For simplicity, we will only look at X — > e + e~, and 
focus on the case that X couples only to electrons and not to protons. To capture a wide 
range of possible "dark boson" scenarios, we allow X to be scalar, pseudoscalar, vector, or 
axial- vector. 4 

In the next section, we summarize the conclusion of our study, that with 1 month to 1 
year of running at the JLab FEL, one can probe an interesting parameter space for the X 
boson. In Sec. 3, we outline our theoretical setup, and review indirect and direct constraints 
on the X boson properties. We study the reach for X in ep collisions in Sec. 4 and show 
how a matrix element method can be used to extend the X boson reach. Comparisons 
to other X boson collider searches appear in Sec. 5. Two benchmark scenarios appear in 
Sec. 6, showing an example analysis strategy as well as a variety of kinematic distributions. 
We consider an alternative displaced vertex search in Sec. 7 and conclude in Sec. 8. 



3 For an alternative low-energy search using a positron beam incident on a hydrogen target, see Ref. [32]. 
4 In the case of a light vector boson, X is often referred to as a U -boson. 
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Figure 1: Summary of this study, showing constraints and reach on the mx-ax plane. The scalar 
couplings arc on the left, vector couplings on the right. The top plots are logarithmic in mx, 
while the bottom plots are linear in mx to highlight the relevant parameter space. Blue curves: 
bounds from lepton anomalous magnetic moments. Orange curves: bounds from prior beam dump 
experiments. Purple curves: luminosity necessary to achieve S/VB = 5 assuming m e + e - invariant 
mass resolution of 1 MeV, and a detector acceptance of —2 < r\ < 2, KE p > 0.5 MeV, and 
KE e ± > 5 MeV. The reach in the pseudoscalar and axial- vector cases are the same as the scalar 
and vector cases, respectively, but the anomalous magnetic moment bounds differ. The points 
labelled "A" and "B" correspond to the benchmark scenarios studied in Sec. 6. 
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2. Summary 

The key results of our study are summarized in this section. We take an electron beam 
with a fiducial energy of 100 MeV, incident on a proton target at rest in the lab frame. 
At the FEL facility, the proton target is likely to be a hydrogen gas storage cell, which 
is sufficiently diffuse to allow recirculation of the FEL beam [34]. In Fig. 1, we show a 
plot of the mx-otx plane for the case that X has scalar or vector couplings to electrons. 5 
Here, mx is the mass of the X boson, and ax = A 2 ^ /Air is the coupling of the X boson to 
electrons, normalized such that it can be roughly compared with «em — 1/137. 
There are three sets of curves shown on the mx~ax plane: 

• Indirect constraints from anomalous magnetic moments. An X boson that couples 
to electrons will radiatively generate a e and, assuming lepton universality, a^. The 
shaded region is excluded, as explained in Sec. 3.1. 

• Direct constraints from beam dump experiments. When the X boson is sufficiently 
long-lived, high intensity beam dump experiments are sensitive to the decay X — > 
e + e~ . The shaded region is excluded, as explained in Sec. 3.2. 

• Discovery reach in ep scattering. As detailed in Sec. 4.2, we use an idealized detector 
with pseudorapidity coverage —2 < r/ < 2 (i.e. tracking up to 15.4° of the beamline), 
kinetic energy thresholds of KE p > 0.5 MeV and KE e ± > 5 MeV, and invariant mass 
resolution of 1 MeV. Assuming no systematic errors, the curves show the integrated 
luminosity needed to achieve a 5a discovery with statistics alone, i.e. S/\H3 = 5 in a 
1 MeV resolution bin centered around mx- If the energy resolution is improved, this 
integrated luminosity required for discovery improves linearly. Also, the reach of the 
experiment can be further improved using a matrix element method, as proposed in 
Sec. 4.3. 

In addition, the points labelled 'A" and "B" indicated the benchmark scenarios considered 
in Sec. 6. 

We see that with an integrated luminosity of 1 ab _1 , there is a range of X boson 
masses and couplings that are consistent with known bounds but visible in low energy ep 
scattering. With an average current of 10 mA, the FEL beam produces 6 x 10 6 electrons per 
second [35], while a hydrogen gas target of thickness 10 19 cm -2 is expected to be technically 
feasible [34]. Thus, the expected luminosity of the JLab FEL setup is 6 x 10 35 cm -2 s -1 , 
which is approximately 1 ab _1 /month. With one month to one year of running, such a 
facility could probe an interesting range of X boson masses and couplings. 

3. New Light Bosons 

We consider four kinds of coupling for the new light boson X: scalar, pseudoscalar, vector, 
and axial-vector. For simplicity, we assume that X only couples to electrons (and other 

The plots for pseudoscalar and axial-vector couplings appear later in the text. The reach for the 
pseudoscalar (axial- vector) is the same as the scalar (vector), but the magnetic moment bounds differ. 
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charged leptons), and show in App. E that proton couplings do not drastically change our 
conclusions. In the scalar and pseudoscalar cases, we augment the standard model with a 
new boson X that couples to the electron field ip e as 

C./p = (As + A p7 5 ) fcX. (3.1) 
For the vector and axial- vector cases, we add a massive vector with couplings 

= + Aa7 M 7 5 ) ^u- (3.2) 

While more exotic operators are possible, Eqs. (3.1) and (3.2) cover the generic possibilities 
for how X can couple to electrons. 

In a complete theory, there is usually some kind of lepton universality, yielding equiva- 
lent couplings of X to muons and taus. Assuming that X does not introduce lepton flavor 
violation, then the coupling of X to the different leptons will either be approximately equal 
or proportional to the lepton masses. 6 As we will see below, though, if the couplings are 
indeed proportional to the lepton masses, then the constraints from a M exclude any of the 
interesting region for this study. 

For convenience where relevant, we define 

X x = ^|A,| 2 + |A p | 2 or vU>| 2 + |A a | 2 , (*x = ^- (3.3) 
In the text, we work in the limit m e <C mx, and present formulas for finite m e in App. A. 
3.1 Indirect Constraints 

Previous studies of the indirect constraints on the X boson appear in Ref. [18]. The 
strongest indirect bounds comes from the effect of the X boson on the anomalous magnetic 
moments of the electron and muon, a e and a,, arising from X boson loops. We take the 
limits to be 5a e < 1.7 • 1CT 11 [19] and <5a M < 2.9 • 1CT 9 [36]. However, it should be noted 
that in the case of vector or pseudoscalar couplings, the bound should not be taken as 
a hard constraint. We have chosen a limit such that the addition of the X boson does not 
significantly change the agreement between experiment and theoretical predictions. Yet 
with the currently observed muon magnetic anomaly the agreement between theory and 
experiment in a M actually improves for vector or pseudoscalar couplings just above the 
constraint. 

If the couplings of the X boson to leptons were proportional to the lepton masses, 
then the constraint would exclude the interesting parameter range for this study. The 
reason is that the production rate for X bosons is proportional to a x , but the constraint 
on a M is on a x = a e x (m^/m e ) 2 . Thus, the effective bound on a e x is almost 5 orders of 
magnitude stronger than if a e x = a x . For this reason, we focus on the case with lepton 
universal couplings. 

6 In the special case that X is a pseudo Nambu-Goldstone boson (such as in Ref. [16]), one expects 
X p = me/ f a , where me is the mass of the lepton and f a is the decay constant. 
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Figure 2: Constraints on the X boson mass mx and coupling ax from anomalous magnetic 
moments, a e = (^j^) e and = We assume lepton universality with a x = a x . 

The formula for the anomalous magnetic moment appears in App. A.l, mirroring 
known results from Refs. [37, 38]. Taking the limit m e <C mx. 

The calculation in the muon case does not admit a simple approximate form since ~ mx 
in the range under consideration, so one must use the full formula from the appendix. 

The constraints on the coupling of the X boson to electrons and muons are shown in 
Fig. 2, assuming a x = a x . For mx > 10 MeV, the bound dominates, giving roughly 
ax 10 -7 - 10 -8 , except for the axial- vector case, for which the bound excludes much 
of the interesting parameter space for this study. Thus, we see that the coupling of the 
X boson to leptons must be 5 to 7 orders of magnitude weaker than the electromagnetic 
coupling «em — 1/137. 

3.2 Direct Constraints 

The strongest direct constraints on the X boson come from beam dump experiments. As 
explained below, other direct constraints from rare meson decay modes, production from 
cosmic rays, or supernova cooling either fall outside the parameter space under considera- 
tion or are subsumed by the beam dump and anomalous magnetic moment constraints. 

Despite the small coupling of the X boson to electrons, it still has a sizable production 
rate in beam dump experiments because of their very high luminosity. As long as the X 
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Beam Dump Constraints (Scalar) 



Beam Dump Constraints (Vector) 




Figure 3: Constraints on the X boson mass mx and coupling ax due to prior beam dump 
experiments, E774 at Fermilab and E141 at SLAC. The red lines correspond to X boson lifetimes. 
The lifetime (and corresponding beam dump constraints) of the pseudoscalar and axial-vector are 
nearly identical to the scalar and vector, respectively, except in the very low mass region. 

boson lifetime is sufficiently long, the decay of X to electrons happens at finite displacement 
from the target. Using the couplings in Eqs. (3.1) and (3.2), the X boson width has been 
calculated in App. A. 2. In the limit m e <C mx, the width is linear with mass 

Roughly, the beam dump constraints are relevant for lifetimes longer than ~ 1CP 3 cm. 

As recently discussed in Ref. [29], the two experiments relevant in the parameter space 
of interest are Fermilab's E774 and SLAC's E141 electron beam dumps, both originally 
searches for MeV-scale axions. E774 consisted of a dump of 5.2 • 10 9 electrons at 275 GeV 
on a 19.6-cm tungsten target, with a 20-cm wide detector 7.25 meters away. With the 
trigger requiring a decay product with energy of 2.75 GeV, 17 events would qualify for 
discovery [28]. The E141 beam dump consisted of a 9 GeV source with 2 • 10 15 electrons 
incident on a 12-cm tungsten target with a detector 35 m away down a 7.5-cm pipe. The 
trigger consisted of a decay dumping 4.5 GeV on the detector. Given the background, 1000 
events would constitute discovery, with the greater number of events necessary due to the 
lack of veto counters on this experiment [27]. 

To determine the excluded parameter space, we mirror the analysis of the vector case 
from Ref. [29] and extend it to the other X boson couplings. In particular, we use an 
approximate formula for the cross section and kinematics for X boson production and decay 
in order to model the coherent nuclear effects involved in the beam dump experiments. The 
excluded regions are shown in Fig. 3, with X boson lifetimes overlayed. The upper diagonal 
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boundary of the excluded region corresponds to when the X boson lifetime is sufficiently 
short that X decays to electrons within the shielding of the beam. Naively, one would 
then expect the upper boundary to be along lines of constant lifetime. The discrepancy 
is due to the fact that lower- mass bosons are created with higher average Lorentz factors, 
allowing more potential X boson events to be seen in the detector downstream past the 
shielding. The lower boundary of the excluded region is determined by the rate of X boson 
production. If the typical X boson decays past the detector, then this boundary would 
be approximately horizontal. The reason is that for an X decay length tx and distance 
to detector L, the fraction of decayed X bosons, approximately L/ix, cancels the l/mx 
mass dependence in the production rate. When the decay length becomes less than the 
distance to the detector, the lower boundary rises at a diagonal. 

Note that the beam dump constraints exclude any of the parameter space accessible in 
this search where X is a very long-lived particle. This is important, because if X were too 
long-lived, then it would decay outside of a typical detector volume, and be inaccessible in 
the ep scattering experiment under consideration. Regions of moderate X lifetime are still 
allowed, so the X might decay with a displaced vertex in ep scattering. We will return to 
this possibility in Sec. 7. 

We can also consider direct constraints that rely on model-dependent assumptions on 
hadronic couplings. Most of these end up being either irrelevant or superfluous in the region 
of interest. For example, B- factories provide constraints on the X boson coupling constant 
though rare T decays, but such constraints currently only apply for nix > 2m^, and are 
thus beyond the range of interest [39] . A recent BaBar analysis [40] does constrain dark 
force models, but relies on a model-dependent signature present only in non-Abelian dark 
sectors. Some proton beam dumps, such as CHARM at CERN [41] do cover some of the 
parameter space in question, but are already subsumed by the electron beam dumps. There 
also exist bounds from the rare pion decay mode 7r° — ► e + e~ [20], which can place some 
additional constraints in a small region close to where the electron and muon anomalous 
magnetic moment bounds intersect. 

There are also potential bounds that end up being irrelevant since the X boson is 
short-lived in the parameter space of interest. For example, X bosons could potentially 
be observed due to cosmic rays interacting with the earth, such that X bosons could be 
seen at detectors such as AMANDA and ICE-CUBE. However, these experiments rule 
out couplings that are smaller (equivalently, lifetimes that are longer) than those we are 
trying to observe here [29] . Additionally, neutrino searches such as LSND and MiniBooNE 
can also be used to place constraints on MeV scale bosons, but again at much lower 
couplings [42]. Another constraint comes from supernovae, where X boson production 
could lead to additional cooling of the core. However, the X boson would have to travel 
at least O(10 km) in order to escape the core. Thus, such constraints require lifetimes for 
the boson to be several orders of magnitude longer than those considered in this search. 
Detailed calculations done for the axion case in Ref. [43] agree with this rough estimate. 

Assuming hadronic couplings, strong bounds on the X boson might be obtained from 
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Figure 4: The envisioned experimental setup to probe e~p — > e~pe + e~ scattering, where a high 
intensity electron beam is incident on a diffuse hydrogen gas target. In order to reconstruct all four 
outgoing fermions, a large tracking volume is needed, and we will consider — 2 < -q < 2 coverage as 
a benchmark. While this is a fixed target experiment where most of the tracks go in the forward 
direction, we take a detector symmetric around r\ = for simplicity. For reference, r\ — 1 is an 
angle to the beamline of = 40.4°, rj = 2 is 6 = 15.4°, rj = 3 is 6 = 5.7°, and 77 = 4 is 6 = 2.1°. 

data mining the extant pion decay data set from the KTeV collaboration. 7 Over the 
course of two runs, approximately 1.8 • 10 6 decays of ir° — ► e + e~7 with invariant mass 
m e + e - > 65 MeV have been reconstructed [44]. Assuming this can be extended to the 
whole invariant mass range, this should give a data set of close to 6 • 10 7 decays. One 
can then look for an X boson by looking for a small peak in the m e + e - spectrum, in a 
search much like the e + e~ — > e + e~7 search described later in Sec. 5. Using an approximate 
formula for the reach from Ref. [21], a search looking for X bosons should be able to detect 
couplings as low as ax ~ few x 10 -8 , an improvement on current bounds, but not matching 
the reach for the search proposed in this paper. 

4. Electron-Proton Scattering Below the Pion Threshold 

We now describe a promising venue for a light X boson search: low-energy ep scattering. 
Despite the above constraints on ax, the cross section for X boson production in ep — > 
ep+X is quite large, on the order of picobarns. The X would then decay promptly to e + e~ , 
yielding an e~pe + e~ final state where the electron/positron pair reconstructs a narrow X 
resonance. The envisioned experimental setup is sketched in Fig. 4, where a high intensity 
electron beam is incident on a hydrogen gas target. Assuming a sufficiently large tracking 
volume, all four outgoing fermions can in principle be reconstructed. 

As we will see, the primary challenge for this search is the large irreducible QED back- 
ground, roughly four orders of magnitude larger than the signal, making this a background 

7 We thank Maxim Pospelov for bringing this possibility to our attention. 



- 10 - 




+ 



Figure 5: Diagrams contributing to the X boson signal. Here, the X boson propagator is evaluated 
in the narrow- width approximation. 

limited search. In this section, we summarize the signal and background calculations and 
present the reach for the X boson, first using simple cuts on the detector geometry and then 
including full information about the event kinematics through a matrix element method. 

4.1 Cross Section Calculations 

In the range of couplings allowed by the anomalous magnetic moment bounds, the width of 
the X boson is an eV or smaller. Therefore, to calculate the signal rate for e~p — > e~p e + e~ , 
we can safely use the narrow width approximation. Since the angular distribution of the 
decay X — > e + e~ is relevant for understanding the vector couplings, we maintain full 
polarization information in the signal process e~p — ► e~pe + e~ as explained in App. B. We 
used a custom matrix element /phase space generator to calculate the signal cross sections, 
and checked the results with CompHEP 4.5.1 [45] . In particular, we used CompHEP to verify 
that any interference between the signal and background processes is a subdominant effect 
given the narrowness of the X boson. 

For the couplings in Eqs. (3.1) and (3.2), the diagrams that contribute to the signal 
cross section appear in Fig. 5. In the limit that m e — ► 0, the cross section for the pseu- 
doscalar (axial-vector) case is identical to the scalar (vector) case. While we keep finite 
m e effects in our calculations, we will only show reaches for the scalar and vector cases, 
since the finite m e effects are small. While there could be contributions to the signal from 
X boson couplings to the proton, we argue in App. E that such effects can be ignored. In 
reconstructing the X resonance, there is combinatoric confusion about which electron to 
pair with the positron, and this confusion is included in our plots. We neglect the electro- 
magnetic form factor of the proton, which is a fair approximation since we are considering 
incoming electron energies E e <C m pro t n- 

The background to e~p — > e~pe + e~ is due to QED radiative processes e~p — > e~p+j* 
with 7* — ► e + e~ and to the Bethe-Heitler trident process, shown in Fig. 6. Details of these 
backgrounds appear in App. C, where we again ignore the proton form factor. We calcu- 
lated the background cross sections using a custom phase space generator interfaced with 
the stand-alone version of MadGraph 4.4. 17 [46], and checked the results using CompHEP. 
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Figure 6: Representative diagrams contributing to the radiative QED background. With indis- 
tinguishable outgoing electrons, there are 12 diagrams in total, consisting of 7* emission off the 
incoming/outgoing electron/proton lines and the Bethe-Hcitler process. In this study, we ignore 
the electromagnetic form factor of the proton, which contributes at most a 5% correction. 

Because we are considering mx ~ E e , thinking of the 7* as coming from initial or final 
state radiation is not a good approximation to the background process. In the region of 
phase space where an e + e~ pair fakes an X resonance, the photon is far off-shell relative 
to the energy scales involved, so we are far away from soft-collinear limit. In particular, 
it should be noted that considerable (constructive) interference increases the background 
above naive expectations from the Weizsacker- Williams approximation. To give a sense of 
how important this interference is, one sees changes on the order of 10% in cross sections 
with e~p collisions versus e~p collisions, whereas the sign of the proton charge would be 
irrelevant in the Weizsacker- Williams picture. 

In principle, there is another background we should consider. Since the proposed JLab 
FEL experiment is really an electron- hydrogen gas collider, one might be concerned about 
backgrounds from e~e~ collisions. In fact, this is only an issue for very low values of mx- 
For an electron beam with energy 100 MeV, the center-of-mass energy of e~e~ collisions is 
around 10 MeV, so radiative M0ller scattering is only relevant for mx < 10 MeV. We saw 
in Sec. 3.2 that such light bosons are already ruled by direct constraints. Moreover, given 
the fact that we imagine using the recoiling proton and electron momentum as a handle 
on the collision process, particle identification should be sufficiently robust to distinguish 
e~e~ from ep collisions. Finally, there are important experimental backgrounds, including 
event pileup and photon conversion, which we do not address in this study. 

4.2 Resonance Reach 

Since the X boson is very narrow, with infinitely good mass resolution and the relatively 
large X boson production rate, one could easily produce the few events required in the 
same e + e~ invariant mass bin to conclude the existence of X. In practice, though, one 
must take into account finite experimental resolution. In Fig. 7, we plot the signal cross 
sections for fixed ax = 10~ 8 as a function of mx, and compare it to the QED background, 
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Cross Sections 




m x (MeV) 



Figure 7: Cross section for the X boson signal with ax — 10~ 8 and the QED background imposing 
that at least one e + e~ pair is in the mx ± 0.5 MeV invariant mass bin. The signal cross section 
scales linearly with ax ■ Generically, the expected signal is four orders of magnitude smaller than 
the background. 

imposing a cut that at least one e + e~ pair within a 1 MeV mass bin around mx- One can 
see that the signal cross section is in the range 10 -2 to 10 2 pb, but the background size 
is generically four orders of magnitude large than signal. Thus, one will need a very large 
integrated luminosity to establish the signal over statistical fluctuations in the background. 

To assess the reach of experiment precisely, one would need to know the true resolution, 
efficiency, and dimensions of the detector. As a rough approximation to the detector 
geometry, we assume full azimuthal coverage, and consider a detector with pseudorapidity 
coverage of 

-2 < r/ < 2, 7/ = -ln^tan0 (4.1) 

(i.e. a tracking system that covers angles as close as 15.4° to the beam line). As we will see 
later in Fig. 8, a more aggressive —3 < i] < 3 coverage (i.e. tracking up to 5.7°) actually 
has a comparable reach. We also impose a constraint on the kinetic energy of the outgoing 
particles: 

KE P > 0.5 MeV, KE e ± > 5 MeV. (4.2) 

While we assume the detector is symmetric about rj = 0, recall that we are considering a 
fixed target geometry, so the tracks dominantly appear in the forward part of the detector. 
Since the QED background has a large forward peak, the pseudorapidity restriction does 
improve the signal to background ratio compared to what is shown in Fig. 7. We only keep 
events where all four outgoing fermions are contained in the tracking volume. 
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Figure 8: Reach plots for variable detector angular acceptance r\. The grey shaded regions cor- 
respond to the indirect and direct constraints from Sec. 3. In all cases, we impose the criteria 
KE p > 0.5 MeV and KE e ± > 5 MeV. We again take an integrated luminosity of 1 ab _1 and 
assume 1 MeV m e + e - resolution. For small values of 77, the reach deteriorates because the signal 
efficiency decreases. For large values of rj, the reach deteriorates because the background has a large 
forward peak. 

As a baseline, we assume that the e + e~ invariant mass resolution is 1 MeV of the 
target mx value. We then calculate the value of ax such that for a given luminosity, 
one can achieve a 5cr discovery with statistical uncertainties alone, meaning we find where 
S/yB = 5 in a 1 MeV mass bin centered on a candidate X mass. Since the background is 
relatively smooth over the kinematic range of interest, the required luminosity for discovery 
scales inversely with the mass resolution: 

C(x MeV resolution) = —C{\ MeV resolution). (4.3) 

x 

The reach for 1 MeV resolution was shown above in Fig. 1, taking 3 MeV < mx < 
100 MeV. 

In Fig. 8, we show how the reach changes as the angular acceptance is varied. Going 
from < 2 to \r]\ < 3 slightly improves the reach for smaller values of mx, though the 
effect is mild. As we will discuss more below, the reason the \r]\ < 2 geometry is so effective 
is that by cutting out the phase space close to the beam line, we decrease the background 
rate without sacrificing much on signal acceptance. 

This S/^/~B procedure to establish reach is only a crude estimate of the true sensitivity. 
In practice, the actual background distribution would have to be fit from the data using 
some kind of sidebanding procedure (see Sec. 6), and one also must pay a trials factor 
in looking for an invariant mass peak since the X boson could be anywhere. With those 
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caveats, we see that with 1 ab 1 of data, one begins to probe the interesting parameter 
regime for the X boson. 

4.3 Reach with Matrix Element Method 

The reach plots in Fig. 1 do not include any kinematic cuts apart from the —2<i]<2 
cut on detector acceptance and the kinetic energy restriction from Eq. (4.2). As we will 
see in Sec. 6, the kinematic distributions for signal events do differ from the background, 
so one might hope that a set of optimized kinematic cuts might improve the reach for the 
X boson. Here we show that a factor of 3 improvement in the reach is in principle possible 
by using complete kinematic information via a matrix element method [47, 48]. 

The matrix element method is often described in terms of a discriminant function [49] , 
but the essential statistics can be understood by considering a weighted measurement. For 
a very narrow resonance X, the signal and background matrix elements for e~p — > e~p e + e~ 
are essentially functions of $4, which is the four-body final state phase space $4 with an 
additional restriction that one of the electron / positron pairs reconstructs a given value of 
mx- For simplicity, we will use the notation $ to refer to $4. 

For a differential signal cross section times luminosity S{§) and differential background 
times luminosity B{&), the naive reach calculation is equivalent to integrating over all of 
<E> with unit weight: 

S = Jd$S($), B = Jd$B{$). (4.4) 

The reach is determined by calculating S/5B = S/yB, where 5 refers to the statistical 
uncertainty in the measurement. Now consider a weighted measurement 

S- eff = J d* B cS = J d$B(<S>)w(<S>), (4.5) 

where w(&) is some weight function. For example, a weight function corresponding to 
hard kinematic cuts is one where u>($) equals either or 1. However, more general weight 
functions still give well-defined measurements. 

The matrix element method calculates the optimal kinematic observable to discrim- 
inate signal from background, which corresponds to choosing the optimal function w(&) 
that maximizes S e ff/6B e g. As derived in App. D, the ideal weighting function is 

= g§f (4-6) 

which yields 

S e fi 



SB cS 



best 



d*|?£. (4.7) 



We can therefore recalculate the reach for the X boson using this ideal value for S e s/5B c s, 
and the results are shown in Fig. 9. As advertised, there is potential factor of 3 improve- 
ment in the reach by using the full kinematic information in the signal and background 
distributions. 
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Figure 9: Reach plots using the matrix element method. The solid curves include the detector 
acceptance cuts —2 < 77 < 2, KE p > 0.5 MeV, and KE e ± > 5 MeV, while the dashed lines have 
no acceptance cuts. The green curves indicate the reach with the matrix element method, and 
the purple curves without. In all cases, we take an integrated luminosity of 1 ab _1 , and the reach 
corresponds to S/V~B = 5 assuming m e + e - resolution of 1 MeV. The matrix element method can 
yield around a factor of 3 improvement in the reach. Note also that the detector geometry alone 
does act like a primitive matrix element method by cutting out the forward region. 

Of course, the matrix element method assumes that the uJbest(^) function is calculated 
using the true signal and background distributions, and this is not possible in practice, due 
to both theoretical uncertainties in the matrix elements and detector effects. Still, one 
might still hope to improve the reach by doing hard kinematic cuts that approximate 
^best(^)- As an example of this, consider Fig. 9. There one sees noticeable improvement 
in S/y/B just from applying the fiducial detector geometry. In this way, the detector 
geometry does act like a primitive In Sec. 6, we will look at lUbest(^) in more detail 

to see what other kinds of hard cuts could be most helpful in teasing out the signal. In 
principle, by using a polarized electron beam, one could obtain additional information from 
the full ep — > ep+X matrix element, but we will not consider polarized beams in this paper. 

5. Comparison to Other Searches 

We argued that low energy ep scattering with at least 1 ab~ x of data was a promising venue 
for looking for a light, weakly coupled X boson. Unlike beam dump experiments where 
the proton recoil spectrum is not measurable, a high intensity electron beam on a diffuse 
gas target allows for full event reconstruction. 

However, ep scattering is certainly not the unique choice of experiment with full re- 
construction potential, and electron-electron scattering or electron-position scattering also 
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have large X boson production rates. Here, we will argue that for the same integrated 
luminosity, ep, e~e~, and e + e~ colliders all offer comparable search power. But given 
the very high instantaneous luminosity achievable at the JLab FEL, we believe that ep 
scattering is favored for X bosons in the range 10 MeV — 100 MeV. 
Consider the following four scattering processes: 

e~p — > e~p + X (fixed target), 

e~e~ — > e~ e~ + X (colliding beams), 

e + e~ — ► e + e~ + X (colliding beams), 

e + e~ — ► 7 + X (colliding beams). (5-1) 

The first one is the ep scattering experiment in this paper, the next two are the equivalent 
processes for e~e~ and e + e~ colliders, and the final search channel is only available for 
e + e~. The distinction between "fixed target" and "colliding beams" is only needed to de- 
termine the relation between the lab frame and the center-of-mass frame of the experiment. 
We assume the same detector technology for all four experiments, with pseudorapidity cov- 
erage in the lab frame of —2 < r/ < 2, the kinetic energy restriction from Eq. (4.2), and 
1 MeV invariant mass resolution. For the e + e~ — ► 7 + X search, one would also need to 
impose a cutoff on photon energy, but as the photons are monochromatic for a given beam 
energy, this would merely correspond to no reach at all once the cutoff is reached and is 
thus not shown. 

In Fig. 10, we show the reach in ax as a function of the available beam energy 

E c s = \fs — m\- 7712, (5.2) 

where y/s is the center-of-mass energy of the collider and rrii are the beam masses. (For ep 
scattering with a 100 MeV electron beam, ~ 95 MeV.) We take a fixed mx = 50 MeV 
and fixed integrated luminosity C = 1 ab _1 . The signal and background cross sections 
were calculated using the same method as in Sec. 4. 

Given the same integrated luminosity and high enough values of E e g, the reach for the 
searches involving four final state fermions are within a factor of 2 of each other. However, 
getting 1 ab" 1 of data in the proposed FEL experiment requires only 1 month of data 
taking, while the maximal luminosity currently achieved in colliding beam experiments is 
1.7 • 10 cm" ~ 2 s^ 1 ~ 0.5 ab -1 /y r - Since high luminosity is critical to probe the parameter 
space of interest, this favors a fixed target experiment for the four fermion final states. For 
the same luminosity, the e + e~ — > j + X search is 1 to 2 orders of magnitude more sensitive 
than the ep search, and we will comment more on this search below. 

It should be noted that existing collider experiments, such as BaBar and Belle, already 
have data sets with integrated luminosities of ~ 1 ab _1 . However, the collisions there 
occur at much higher energies. At those energies and for the same search strategy, many 
additional backgrounds, both reducible and irreducible, are present, increasing the difficulty 
of the analysis. In fact, studies of these detectors typically focus on the case with mx > 
2m as the decay of the X boson to muons is more easily reconstructed than the decay to 
electrons [24, 39]. 
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Figure 10: Reach for an 50 MeV X boson, varying E c g — y/s — rni — m-x. We compare the present 
ep scattering proposal with three alternative searches: e~e~ — > e~e~ + X , e + e~ — > e + e~ + V, and 
e + e~ — > 7 + V. In all cases, we assume a detector acceptance of —2 < 77 < 2, KE p > 0.5 MeV, and 
KE e ± > 5 MeV. While the in principle reach in e + e~ — * 7 + X appears to be better than the other 
three, we discuss the challenges of that search in the text. For the same integrated luminosity and 
large enough E c g, the three searches with four outgoing fcrmions have comparable reach, which 
favors the ep scattering proposal where high luminosity is more readily achievable. 



On the other hand, even at the higher energies of existing data sets, the e + e~ — ► 7 + X 
channel does look quite promising for an X boson search due to the lack of hadronic back- 
grounds. That said, two factors lead this channel to be more complimentary than compet- 
itive with the search we propose. Our naive estimate of the 7 + X reach assumed perfect 
reconstruction of every event meeting the detector geometry cuts. This is significantly 
more difficult with a search using on-shell photons as there is no longer tracking informa- 
tion for every particle. In order to identify the single energetic photon in the event with 
high accuracy, it is necessary to put tighter cuts on the photon angle to get farther away 
from the beam pipe. Estimates in Ref. [22] indicate that when such cuts are put in place, 
the reach in coupling actually becomes comparable to that of the ep search. Additionally, 
the search for 7 + X is complicated by photon conversion in the tracking volume from the 
much larger e + e~ — > 77 process. While this can be offset by cutting on displaced vertices, 
such an approach becomes difficult for e + e~ mass bins below 50 MeV, and impossible in 
mass bins below 20 MeV [50]. For X bosons above 50 MeV, one should be able to extend 
the 7 + X search already done for X — > [40] to X — > e + e~, making the search 

complementary to ep scattering by filling in the mass range between ~ 50 MeV and 2m M . 

The preceding discussion of the benefits of ep scattering applies for X bosons in the 
mass range 10 MeV < mx < 100 MeV. For higher values oimx, the JLab FEL simply does 
not have the kinematic reach achievable in other colliders. The maximum sustainable FEL 



- 18 - 



Reach for 1 ab 1 (Scalar) Reach for 1 ab 1 (Vector) 




m x (MeV) m x (MeV) 



Figure 11: Reach plots for ep collisions, increasing the electron beam energy from E e — 100 MeV 
to E e — 140 MeV, the maximal sustainable beam energy available at the JLab FEL. Assumptions 
about detector geometry, integrated luminosity, and energy resolution are the same as the previous 
figures. 

beam energy is around 140 MeV [34], limiting the in-principle reach to mx < 131 MeV. 
Moreover, even if one were able to get higher energy electron beams, it is no longer clear 
whether ep scattering would pose any advantage, since for electron beam energies above 
the pion mass, inelastic scattering channels open up, increasing the number of tracks per 
beam crossing. 

Below the pion mass threshold, though, the reach in ep collisions can be improved 
in going to somewhat higher E e Q as seen in Fig. 10. A more detailed look is shown in 
Fig. 11, which compares the reach for E e = 100 MeV to the maximum sustainable FEL 
energy of E e = 140 MeV. For X boson masses close to the kinematic limit, the higher 
energy electron beam gives improved sensitivity, though at low masses, most of which have 
already been ruled out by beam dump constraints, the reach gets slightly worse. 

6. Benchmark Studies 

In this section, we consider two benchmark X boson scenarios, indicated on Fig. 1, 

A : m x = 50 MeV a x = 10" 8 , (6.1) 
B : m x = 20 MeV a x = 3 • 10" 9 , (6.2) 

with both scalar and vector couplings in each case. These points were chosen to be roughly 
close to the 1 ab -1 reach lines in Fig. 1. We will show an example analysis strategy for X 
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Figure 12: Invariant mass distribution for m e + e - , comparing the QED background to benchmark 
models A (left) and B (right). The fluctuations in the background distribution at low m e + e - are 
from Monte Carlo statistics, and are not indicative of expected statistical fluctuations. The signal 
includes the combinatoric background from pairing the "wrong" electron with the positron. These 
plots include the detector acceptance criteria —2 < r\ < 2, KE p > 0.5 MeV, and KE e ± > 5 MeV. 
Note the four orders of magnitude difference between the expected signal and background. 

boson signal extraction and then show various kinematic distributions to highlight which 
parts of phase space are most sensitive to X boson production. 

To begin, consider reconstruction of the X boson resonance. In Fig. 12, we show the 
invariant mass distribution for e + e~ pairs, taking the fiducial detector acceptance from 
Eqs. (4.1) and (4.2). Since the final state electrons are indistinguishable, the plot includes 
two histogram entries per event. In the case of the signal distribution, there is the expected 
spike at mx accompanied by a combinatoric background from pairing the positron with 
the "wrong" electron. We see that there is a four orders of magnitude difference between 
signal and background, consistent with Fig. 7 which did not include any detector acceptance 
effects. 

To show how the X boson resonance could be seen despite the large background, 
Fig. 13 shows a simulated distribution of m e + e - created as follows. First, we generate 1 
ab" 1 of background pseudo-data for the m e + e - distribution and add it to 1 ab _1 of signal 
pseudo-data. We take the combined signal plus background distribution, and fit it to an 
ad hoc functional form: 

= AT(m) a (m max - m)\e m ) c , (6.3) 

dm 

where N, a, b, and c are fit coefficients, and m max is the maximum kinematically allowed 
value for m e + e ~. We then plot the fractional difference between the pseudo-data and the 
final fit function in Fig. 13. For the case of the vector couplings where the benchmark points 
lie above the 1 ab _1 reach line, a peak at the m e + e - distribution at mx is indeed visible 
above the expected statistical fluctuations in the background, showing that a sideband 
procedure for extracting the background is feasible. For the scalar benchmarks, no such 
peak is visible, as expected since these benchmarks lie below the 1 ab _1 reach line. 
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Figure 13: Simulated m e + e - distributions compared to a background fit. This plot was made with 
1 ab -1 of signal and background pseudo-data, assuming 1 MeV invariant mass resolution. The blue 
curves show the expected statistical uncertainties in the background, la (solid) and 5a (dashed). 
For the vector couplings (green), both model A (left) and B (right) show a prominent bump in the 
dilcpton invariant mass distribution, which is expected since these benchmark points lie above the 
1 ab -1 reach line. For the scalar couplings (red), a bump cannot be seen, since these benchmark 
points lie below the 1 ab -1 reach line. 



The X boson resonance signal could be further enhanced by using the matrix element 
method of Sec. 4.3, which would be applied at each value of m e + e - . To see how the matrix 
element method would affect our two benchmark points, we plot some example kinematic 
distributions in Figs. 14, 15, 16, and 17. Except as indicated, these plots are made using 
the fiducial detector geometry cuts of — 2 < rj < 2, KE P > 0.5 MeV, and KE e ± > 5 MeV, 
and require that at least one electron/positron pair reconstructs mx- The electron that 
reconstructs the candidate resonance is the active electron, and the other is the spectator 
electron. In addition to the raw background and signal cross sections, we also plot the ideal 
weighting function from Eq. (4. 6). 8 Large values of the weighting function correspond to 
regions of phase space that have the highest sensitivity to the X boson. 

In Figs. 14 and 16, we show the momentum and angular distributions for the four 
outgoing fermions. All of the detector geometry cuts are imposed, except that the cut 
corresponding to the plotted distribution is indicated via shading. 

From the momentum distributions, we see that the kinetic energy cuts from Eq. (4.2) 
do not cut out much of the most sensitive region, and indeed the cut on proton and 
spectator electron kinetic energy enhances the signal relative to the background. For the 
angular distributions, because the background has a strong peak when an electron scatters 
in the forward direction, the weighting function is suppressed near cos(f? c i cctr0 n) = 1- This 
explains why in Fig. 8, the reach did not improve much in going from an \r/\ < 2 cut to an 
1 77 1 < 3 cut. Away from this peak, the weighting function has roughly flat sensitivity to the 

8 As a sanity check, we verified that the weighting function is approximately equal to the ratio of the 
signal and background in a given observable bin, appropriately normalized to the phase space volume of 
that bin. See App. D. 
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electron/positron angle, indicating that large angular acceptance is important for X boson 
reconstruction. Because of energy-momentum conservation, the proton can only scatter in 
the forward direction, and the most sensitive region is in fact in the most forward region. 

In Figs. 15 and 17, we show the pair-wise invariant mass distributions, as well as three 
angular distributions. The out-of-plane angle is the angle the reconstructed X boson moves 
relative to the plane defined by the incoming electron beam and the spectator electron. 
The decay angle is the angle between the X boson momentum and the outgoing positron 
momentum, as measured in the X boson rest frame. Also plotted is the angle between the 
spectator electron and the reconstructed X boson. 

The invariant mass distributions do show some ability to distinguish signal from back- 
ground. Especially promising is the invariant mass between the spectator electron and the 
positron, which peaks at small values for the background. The angular distributions show 
less promise, as the ideal weighting functions are relatively flat. Note that the weighting 
function for the scalar and vector cases do have different shapes in the angular distribu- 
tions, which explains why the reach plots in Sec. 4 have different mx dependence even 
after correcting for the total signal cross section. 

Once the X boson is discovered, the decay angle could be useful for distinguishing the 
scalar coupling from the vector case. This distribution is flat for the scalar signal, but 
encodes non-trivial angular information in the vector case. This angular variation is much 
smaller than the background, though, so a careful analysis would be necessary to extract 
the nature of the X boson coupling. 
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Figure 14: Momentum and angular distributions for model A, with rax = 50 MeV and ax = 10 -8 . 
The QED background is restricted to have one e + e~ pair reconstruct rax, and the corresponding 
electron is called the active electron while the other is the spectator electron. These plots include 
detector acceptance cuts, but the cut corresponding to the plotted distribution is indicated by 
shading. The solid blue curves are the QED background, and the solid red (green) curves are 
the scalar (vector) signal. The dashed red (green) curves are the ideal weighting functions for the 
scalar (vector) case with arbitrary normalization, which are large in the region of phase space most 
sensitive to X boson production. 
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Figure 15: Additional kinematic distributions for model A, with the same criteria and labeling as 
Fig. 14. Shown are five pairwise invariant mass distributions (the active electron-positron invariant 
mass would of course just give a peak at n%x). The out-of-plane angle is the between the recon- 
structed X boson and the incoming electron/spectator electron plane. The decay angle is the angle 
of the X boson decay products relative to the X boson momentum, measured in the X rest frame. 
Also shown is the angle between the spectator electron and the reconstructed X boson. 
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Figure 16: Momentum and angular distributions for model B, with mx = 20 MeV and ax = 
3 • 1CP 9 , analogous to Fig. 14. 
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Figure 17: Additional kinematic distributions for model B, analogous to Fig. 15. 



- 26 - 



Vertex Reach for 1 ab 1 (Scalar) Vertex Reach for 1 ab 1 (Vector) 




m x (MeV) m x (MeV) 



Figure 18: Reach in ep scattering using a displaced vertex search strategy. The curves correspond 
to producing 10 signal events with an integrated luminosity of 1 ab -1 . The purple curves assume 
that the total displacement between the interaction vertex and the X boson decay vertex can be 
measured. The red curves assume that only the transverse displacement from the beam axis can be 
measured. Note that the plotted range for ax differs from the previous reach plots in this paper, 
and that an additional bound from the E137 experiment is shown. 

7. Prospects for Displaced Vertices 

We saw in Sec. 3.2 that there is a region of parameter space where the X boson has not yet 
been ruled out by the beam dump experiments but has a lifetime long enough to leave a 
displaced vertex in ep collisions. While the X bosons produced in low energy ep collisions 
are generically not very boosted, there is enough of a tail in the lifetime distribution that 
a few X bosons will have large displacements. Moreover, with a diffuse gas target, one can 
in principle probe smaller vertex displacements than in a solid target. 

Since the QED background does not generically lead to displaced vertices, this suggests 
that the reach could be extended beyond that presented in Sec. 4. Of course, there are 
instrumental effects that can lead to fake displacement, but we ignore these in this study. 
For the signal, the displacement can be straightforwardly calculated by convolving the X 
boson momentum spectrum with the exponentially falling lifetime curve, properly taking 
into account boost factors. 

With no expected irreducible backgrounds, we estimate the naive reach by considering 
the region of parameter space where at least 10 displaced vertices would be observed. 
The expected reach using two different search strategies is presented in Fig. 18. The 
more aggressive strategy attempts to reconstruct the total displacement of the X boson. 
That is, the detected recoiling proton and spectator electron define the interaction vertex, 
and electron-positron pair from the X boson define the decay vertex. We plot the reach 
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assuming a 1 mm or 0.1 mm total displacement could be observed. Since there can be 
multiple scatterings per beam crossing, this method faces a background from uncorrelated 
scattering events. 

An alternative strategy is to merely reconstruct the transverse displacement of the 
decay vertex from the beam axis. Since the X boson momentum is peaked in the forward 
direction, this noticeably reduces the reach at the same displacement resolution. Assum- 
ing transverse vertex displacement is easier to detect than total displacement, one might 
imagine that increased resolution might compensate to make this the preferred strategy. 
Note that the transverse displacement search faces a possible background from photon 
conversion within the gas target. 

Both displaced vertex strategies probe a different part of the parameter space from the 
search in Sec. 4, and therefore can be viewed as complementary to the QED background- 
limited analysis. In particular, X bosons with couplings two order of magnitude smaller 
than those accessible with a direct bump hunt could be seen. At such small couplings, other 
prior beam dump experiments come into play beyond those discussed earlier in Sec. 3.2. In 
Fig. 18, the beam dump constraint at the lowest couplings comes from the E137 experiment 
at SLAC [51], calculated analogously to the constraints in Sec. 3.2. 

8. Conclusions 

Low energy electron-proton scattering is one of the basic processes in the standard model. 
It is therefore intriguing that new physics might be discovered in a regime that is thought to 
be dominated by elastic and quasi-elastic QED processes. Recent astrophysical anomalies 
have motivated a new paradigm for dark matter, where heavy dark matter interacts with 
a light, weakly coupled boson. New low energy, high intensity scattering experiments are 
an ideal setting to constrain (or confirm) this exciting scenario. 

We have argued an X boson with ax ~ 10~ 8 and 10 MeV < mx < 100 MeV could 
be discovered in low energy ep scattering with around 1 ab" 1 of data assuming 1 MeV 
invariant mass resolution. Since the search for X — > e + e~ is background limited, it is 
crucial to have an experiment with good energy resolution and very high statistics. We 
believe that the unique combination of high luminosity with full event reconstruction makes 
this a compelling experimental proposal for the JLab FEL. This proposal is complementary 
to the beam dump experiments envisioned in Refs. [21, 29], which are better suited for 
smaller value of ax and larger values of mx- 

We have shown that a matrix element method which uses complete kinematic infor- 
mation about the signal and background can increase the sensitivity to the X boson by 
about a factor of 3. Though not studied in this paper, a polarized electron beam could be 
useful in extracting additional matrix element information. Of course, the most straight- 
forward way to increase the sensitivity of the experiment is to improve the invariant mass 
resolution beyond our fiducial value of 1 MeV. Finally, if reconstruction of 1 mm or 0.1 
mm displaced vertices are possible, then ep scattering could probe an interesting region of 
X boson parameter space with smaller couplings. 
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A. Finite Mass Calculations 
A.l Anomalous Magnetic Moment 

A new light boson, whether it is a scalar or vector, will contribute to the anomalous 
magnetic moment (a = of leptons at the one- loop level. In the low masses being 

considered for the X boson, limits on new contributions to the moment are the main 
indirect constraint and must be calculated accurately. In the scalar/pseudoscalar case, a 
convenient choice of parameterization yields [37] 



Sa '/p = 77^2 / dz ~ ~2 , n Ivrb (A- 1 ) 
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This final form only holds in the limit mg <C mx, and is thus not appropriate for o„. 

For the case of a vector coupling, the term proportional to p^p u in the numerator of 
the propagator drops out of the calculation entirely due to the Ward identity. However, for 
axial coupling this is not the case. At the same time, the extra factor of \jm\ this term 
introduces means it only contributes at 0(m^/m x ). Thus, it can be safely ignored in the 
electron case, but for constraints from the muon anomalous magnetic moment it must be 
included for accurate results. The full calculation yields [38] 
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where again the approximate form only holds for mi <C mx- 



A.2 Lifetime 

If mx < 2771^, then for the couplings envisioned in Sec. 3, the only allowed decay mode 
of the X boson is X — > e + e~. In more general scenarios, the X boson might decay to 
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Figure 19: Curves of constant X boson lifetime on the ax vs. rax plane. 

neutrinos or photons or other new light fields, but such final states face very different 
backgrounds and are outside the scope of this study. 

The tree-level widths from the couplings described in Sec. 3 are 



r*/ P = ^[(K + K) - K^r ) \l 1 - 7-r> ( A - 5 ) 




r u / a = ^ ( (Xl + Xl) + 2{Xl - 2\l)^t) \h ~ % (A.6) 
' Hit \ mir I \ mi 



x 



In the limit m e < mj, these simplify to T s / p = ^(A^ + A^) and T v / a = ^(A^ + A^). A 
plot of the X lifetime is shown in Fig. 19 with full m e dependence, which introduces small 
differences in the scalar /pseudoscalar and vector /axial- vector cases at masses close to the 
electron mass. 



B. Signal Calculation 

In this appendix, we review the narrow width approximation and apply it to X boson 
production. The signal for X boson production is e~p — > e~p + X — > e~pe + e~ . While 
other diagrams with internal X boson propagators can contribute to the e~pe~e + final 
state, only the diagrams shown in Fig. 5 yield a narrow peak in the e + e~ invariant mass 
distribution near mx- 

It is well-known that processes involving a resonance with a small width can be treated 
by the so-called narrow width approximation (see Ref. [52] and references therein). In the 
case of X boson production, the idea is to write the square of the full two-to-four matrix 
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element of diagrams in Fig. 5 as 

|M(2 - 4)|* = | *!>(,, ...) x D (,), D(4 ) = tf _ m ^ + r>r (B.1) 

where g is the four-momentum flowing through the X boson propagator. One then ap- 
proximates the square of the X boson propagator as 

D{q)~^— 5(q 2 -m x ). (B.2) 
m x T x 

Using the cluster decomposition of four-body phase space into a product of three-body and 
two-body phase space, 

d$ 4 = id$ 3 dg 2 d$ 2 , (B.3) 
the fully differential cross section takes the form 

Signal = J j d$ 4 ^|M(2^4)| 2 



spins 



1 1 



/^ 3 /d<i> 2 ^|Af(g,...)| 2 ^. (B.4) 



spins 



Here, F is the incoming flux, which in the case of a fixed target experiment with incoming 
electron energy E e equals 

F = AE e m p . (B.5) 

When the X boson has scalar or pseudoscalar couplings, the matrix element of the full 
two-to-four process takes a convenient factorized form: 

\M(q, ...) s \ 2 = \M S (2 3)| 2 \M S (1 2)| 2 . (B.6) 

For calculating the total cross section, one can perform the d3>2 integral in Eq. (B.4) 
analytically: 

Signal scalar = jj d* 3 | ]T \M S (2 - 3)| 2 (^g— ) ■ M 

spins 

The factor in parentheses is just the branching fraction of X to the e + e~ final state. Of 
course, for any real observable, there are always cuts present, at minimum on the detector 
geometry. So in practice, one must use the full formula in Eq. (B.4) and multiply the 
integrand by the desired observable function. 

When the X boson has vector or axial- vector couplings, the matrix element does not 
factorize, but the amplitude does factorize into a contraction between the production of 
the resonance, the decay of the resonance, and the numerator of the X boson propagator: 

\M v (q, ...)| 2 = |M£(2 -> 3)M£(1 - 2) (g, u - |j) | 2 . (B.8) 
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It is worth mentioning that in the vector case, the q^qu term in the propagator vanishes 
because of Ward identity: 



since the external electron/positron are on-shell. In the axial- vector case, the Ward identity 
is no longer true, but because of chiral symmetry, the correction from the q^q v term is 
suppressed by powers of m e and therefore small. Because of the X boson propagator 
factor, there are non-trivial angular correlations in the X boson decay. 

There is one subtlety in using the narrow width approximation with identical particles 
in the final state. Experimentally, we cannot determine which final state electron came from 
the decay of the resonance and which one is the scattered incoming electron. Therefore, in 
principle, we must add to the diagrams in Fig. 5 another two diagrams with the electron 
legs in the final state interchanged. However, the effect of adding these diagrams but 
including a symmetry factor of 1/2 into the phase space for identical electrons leaves 
Eq. (B.4) unchanged. Since Eq. (B.4) is valid for all regions of phase space, the differential 
distribution in m e + e - will include not just a delta function spike at mx, but also the correct 
combinatoric background. For any real experimental observable, this delta function spike 
will be properly smeared out by the experimental resolution, as long as the experimental 
resolution is coarser than Tx- 

C. Background Calculation 

The background to X boson production consists of quasi-elastic QED interactions, with 
example diagrams given in Fig. 6. Since we are considering incoming beam energies below 
the pion mass, no QCD interaction are relevant, and the proton remains intact after being 
struck by the electron. This simplifies the analysis of the standard model background to 
the proposed signal dramatically. Also, since our target is hydrogen gas, we do not have 
to worry about nuclear excitations. 

Given that the background is QED, we can safely consider tree-level diagrams alone, 
since loops are suppressed by the small electromagnetic coupling. That said, while the 
one-loop corrections to the background are much smaller than the tree-level result, they 
are also expected to be much larger than our signal. However, one-loop corrections are 
not expected to dramatically change the shape of the background (and certainly not give a 
peak at finite m e + e -), so for the purposes of extracting a signal peak in the invariant mass 
distribution, the tree-level background result will suffice. 

The largest correction we are neglecting comes from the electric form factor of the 
proton 



which had been included in the study in Ref. [30]. With mo ~ 700 MeV and q < E e = 
100 MeV, this will at most yield a 5% change in the background calculation. Like the 
one-loop effects, though, the electromagnetic form factor is not expected to dramatically 
change the background shape. 



q v M^{l -> 2) = 





(C.l) 
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There are twelve QED diagrams that contribute to the process e~p — > e~pe + e~ . Just 
as in M0ller scattering, one must be mindful that the interchange of identical fermionic legs 
adds an additional minus sign to the Feynman rules [53]. The background cross section 
equals 

<7background(0) = Jp ^ 4 24 5^ I ^background 1 2 C(^4) , (C2) 

spins 

where 1/2 is a symmetry factor from having two identical electrons in the final state, 1 /4 
the average over initial polarizations and 0(^4) is an arbitrary observable. 

In particular, for calculating the background in the signal bin for a given value of mx, 
we use the theta function 

0($ 4 ) = 9 {{m x - Am/2) 2 < q 2 < (m x + Am/2) 2 ) , (C.3) 

where q 2 is the invariant mass of an outgoing e + e~ pair, and Am is the invariant mass 
resolution. Because the background is a steeply falling distribution in q 2 , as shown in 
Fig. 12, this way of calculating the binned background gives slightly more realistic values 
than 

0{<S> 4 ) =2m x Am5{q 2 -m 2 x ), (C.4) 
though both measurements agree in the small Am limit. 

D. Matrix Element Method 

In this appendix, we derive the ideal weighting function to be used in the matrix element 
method from Sec. 4.3. To start, consider an unweighted measurement 

S = Jd$S(<5>), B = Jd$B(<5>). (D.l) 

The statistical uncertainty in the background 5B can be determined in terms of the Poisson 
uncertainty at each point in phase space 5B(&), 



8B($) = y/B($), 5B = J J d$ [6B(<5>)] 2 = ^/B, (D.2) 



and we recover the familiar formula that S/5B = S/\^B. 9 

Now imagine doing a weighted measurement over phase space: 

S eS = /d$5($)w($), B cS = J d<S>B($)w($). (D.3) 

The statistical uncertainty in B e s is 



SB cB = J j d$ [5B(<S>)} 2 w(<S>) 2 = J J d$B(<S>)iu($) 2 . (D.4) 



9 Strictly speaking, we should really consider statistical uncertainties in both the signal and background, 
but the background is so much larger than the signal that this is superfluous. 
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Note that the reach in S e s/5B e Q is independent of the normalization of w($>). 

To find the ideal measurement function, we simply need to use a variational method 
to solve for the condition 

d / See \ 



dw($>) \5B c{ { J 

Using the fact that 



0. (D.5) 



d d B($)w($) 

SeS = S(^), -^—5B eS = > , (D.6) 



dw($) v " dw($) 5B eS 

it is straightforward to solve Eq. (D.5) to find 

Whes * m ~ s cS W) bWY ( } 

where we have used the fact that the overall normalization of lUbcstC^*) ^ s irrelevant for 
determining the reach. Plugging Wbest(^) into S e g/5B e s, we find 



SB^s 



best 



d*S££. (d. 



For S($>) proportional to B($) (i.e. no kinematic shape differences between the signal and 
background), this formula reverts to the standard S/yB. 

When we plot the weighting function in Figs. 14, 15, 16, and 17, we are actually 
plotting 



/ 



d<S>w hcst ($)(D(<S>), (D.9) 



where 0($) is the observable corresponding to a histogram bin. As a cross check of the 
weighting function, we checked that for most observables, this function is well-approximated 
by 

/d$5($)0($) 



j d$0($), (D.10) 



Jd$ £($)£>($) 

i.e. the binned signal over background ratio corrected by a phase space volume factor. 



E. Generalized Couplings 

In order to study the X boson phenomenology in a model-independent way, we assumed in 
Sec. 3 that the X boson only had couplings to electrons and not to protons. Here, we relax 
this assumption within the context of several proposed models, to see how the X boson 
reach using ep scattering is affected. We find that, for varying reasons, couplings to the 
proton can be ignored. 
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E.l Kinematic Mixing 

In some of the best motivated dark force scenarios, the couplings of the X boson are 
proportional to the electromagnetic couplings [14, 15]. This occurs when the X is a vector 
boson that kinematically mixes with the photon. That is, the vector X^ couples directly 
to the electromagnetic current, albeit with a suppression factor e, 

This yields ax = e 2 OEMi where in concrete models, typical e values are 10 -3 to 10~ 4 . 10 

Adding the X boson coupling to the proton allows for additional diagrams for X boson 
production. Therefore, there is the potential for noticeable contributions to the signal cross 
section from the proton coupling, whether through a direct contribution to the cross section 
or through constructive or destructive interference. By explicit computation, however, we 
have checked that the proton coupling is largely irrelevant, and the essential physics can 
be understood by working out the kinematics of the relevant situations. 

The two differences in producing an X off the photon or electron are the presence of an 
electron versus a proton propagator and the momentum transfer through the exchanged t- 
channel photon. If one notes that neither the proton nor the X boson are very relativistic in 
the energy range under consideration, one can approximately say that E p ~ m p , Ex — mx, 
and E e ~ mx at any point along the electron line. Then, the fermion propagator can be 
shown to be 0(1/ mx) in both cases. However, by seeing how momentum flows through 
the i-channel photon propagator, one can show that in the case of X boson production off 
the electron line, \t\ ~ 0(m 2 x ), while off the proton line, |i| ~ 0(mxm p ). Thus, the cross 
section for X boson production off the proton line is suppressed by 0(m x /m p ). 

This suppression can be understood intuitively as arising from the fact that with a 
proton at rest, enough energy needs to be exchanged in the t-channel to actually create an 
X boson, while with the electron already having sufficient energy, one only needs enough 
momentum transfer to move the electron into a kinematically valid region for X radiation. 

One might worry that because of interference terms, the suppression would only be 
0(mx/m p ). However, in addition to the mass dependence, the momentum exchanged 
through the t-channel photon has angular dependence. This is minimized when the X 
boson is produced collinearly with the fermion line off which it is produced. This means 
that the matrix elements for X production off the electron and off the proton are peaked 
in entirely different regions of phase space, producing little overlap. The interference terms 
end up being suppressed by another 2 orders of magnitude when the diagrams are explicitly 
calculated and integrated over. Thus for mx < 100 MeV, the corrections to the cross 
sections from including proton couplings are around 1%, and do not measurably change 
the reach plots displayed earlier. 

E.2 Axion-Like Coupling 

An alternative framework for the X boson is where the dark sector couples to the standard 
model through a pseudo Nambu-Goldstone boson, termed an axion portal [16]. In that case, 

10 In these models, the typical X boson mass is around 1 GeV. However, models with lighter X bosons 
closer to 100 MeV are still plausible. 
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one expects a coupling of X p = mi/ f a for elementary fermions, where f a is the axion decay 
constant. In such a setup, the constraints from the muon anomalous magnetic moment 
completely rules out the region of electron couplings that might be probed in our setup, as 
discussed in Sec. 3.1. 

However, most axion portal models predict couplings of X p ~ mj/f a for composite 
particles like the proton. Thus, with the coupling to the proton O(10 3 ) greater than to the 
electron, one might think it possible to still see a signal from X boson production off the 
proton, despite the kinematic suppression discussed above. It turns out, though, that a 
decay constant small enough for this to be possible has already been ruled out by K decay 
branching ratios [54], which placed a lower bound of f a = 0(100 TeV) for axion masses 
< 2\i. This leads to a coupling to the proton of at most ax = O(10~ 10 ) which, combined 
with the kinematic suppression, would not be detectable with this search. 
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